Quantum revivals and carpets in some exactly solvable systems 
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We consider the revival properties of quantum systems with an eigenspectrum En oc n^ , and 
compare them with the simplest member of this class - the infinite square well. In addition to 
having perfect revivals at integer multiples of the revival time ta, these systems all enjoy perfect 
fractional revivals at quarterly intervals of t_R. A closer examination of the quantum evolution is 
performed for the Poschel- Teller and Rosen-Morse potentials, and comparison is made with the 
infinite square well using quantum carpets. 
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Over the past ten years or so there has been a growing 
interest in the quantum dynamics of simple systems, mo- 
tivated in part by the richness of new phenomena such as 
revivals ||l| and quantum carpets Ig-Ql. Indeed, the phe- 
nomenon of revivals is not merely a theoretical construct, 
but has been observed in ion traps Q, Rydberg atoms 
PI , and semiconductor wells m ; and has been used to dif- 
ferentiate ionization pathways in potassium dimers [Bj. 
Quantum revivals are similar, but distinct from, quan- 
tum Poincare recurrences, which have been studied re- 
cently in the context of the kicked rotator |H|. In the 
former, one is interested in the deterministic reconstruc- 
tion of the wave-function during its evolution inside a 
fixed potential; whilst in the latter, interest is focused on 
the decay of the return probability in "mixed" regions 
of phase space, as a measure of the quantum chaos in 
the (usually forced) system. For the most part, analytic 
studies of revivals and carpets have concentrated on the 
infinite square well (ISW) potential, which is known to 
have perfect revivals and fractional revivals |1C[| , and also 
a quantum carpet composed of rays (straight lines in the 
space-time plane) [^ . These properties have been under- 
stood on the basis of the quadratic dependence of the 
energy on quantum number, along with the fact that the 
eigenf unctions are elementary trigonometric functions. 

It is well understood that perfect revivals can only 
occur for systems whose energy spectrum is purely 
quadratic in the quantum number R] . If the dependence 
is purely linear (harmonic oscillator) the only time scale 
is the classical period of oscillation, while for more com- 
plicated energy spectra, revivals will be imperfect due to 
modulations from the super-revival time scale. Although 
many studies have been devoted to the simplest quan- 
tum system with a quadratic energy spectrum - namely, 
the ISW - there has been less attention paid to the host 
of other potentials which share this property (although 
see Ref. [O for a discussion of the autocorrelation func- 
tion for the Morse potential). It is guaranteed that these 
systems will have perfect revivals, but what can one say 
about fractional revivals, and the existence of quantum 
carpets {i.e. hidden structures in the space-time plot of 
the probability density)? 



We shall begin with some very general remarks about 
fractional revivals. Consider a system with a purely 
quadratic, nondegenerate energy spectrum En = a^v?, 
(n = 0, 1, 2, • • ■), and with a potential V{x) centered at 
X — Q. We take the potential to be an even function of 
X. In this case the eigenfunctions (j)n{x) will have a def- 
inite even or odd symmetry, alternating as the quantum 
number increases, the ground state naturally being even, 
since it has no node. So we have 0„(— x) = (— l)"(/)„(a:). 
We prepare the wave function of the system in an initial 
state specified by the energy eigenfunction expansion 



1p{x,0) = y^^Cn(f>n{x) 



(1) 



We restrict ourselves to contributions from bound states 
only. The time-evolved wave function is given by 



ip{x,t) = ^c„(/)„(a;)exp[-i£'„t] , 



(2) 



where we have chosen units of ?i = 1 . Given the quadratic 
dependence of the energy levels on n, it is easy to see 
from Eq.(g) that the wave function will be identical to 
its initial state at integer multiples of the revival time 
tR = 2TT/a\ 

Now consider the wave function at a time equal to one 
half of tu. One easily finds 



V'(a;,tfl/2) = ^c„(/i„(a;)exp[-i7m^] . 



(3) 



Given that e *^" == (—1)^, we have 

V'(a;,tfl/2)= ^ c„(/)„{x) - ^ Cn(f>n{x) . (4) 



Ddd 



Returning to the initial wave function, one may use the 
parity properties of the eigenstates to demonstrate that 



■0(0;, 0)= ^ Cn(j>n{x) + ^ C„^„(x) , 
n even n odd 

iIj{-x,0) ^ ^ Cn(f>n{x) - ^ c„0„(a;). (5) 



3dd 



Clearly, on comparing Eqs. (Q) and (^) we have the per- 
fect fractional revival ijj{x,tR/2) = ijj{—x,0). This result 
may appear to follow from the symmetry of the potential 
and time reversal invariance; however, this is not the case 
[cf. the discussion following Eq.(|l^)]. 

A less obvious result follows, however, when we study 
the wave function at one quarter of the revival time. We 
have 



V'(a;,ifi/4) 



Z^ 



Cn4>n{x) exp[—inn^/2] 



(6) 



Considering the phase for n = 0, 1,2,3 mod(4) one can 
easily establish that 

ip{x,tR/A)^ ^ Cn4>„{x)-i ^ c„0„(x) . (7) 

n even n odd 

Solving the two expressions in Eq.(||) for the odd and 
even sets of modes, we find the perfect fractional revival 

V(x, tn/i) = ^^i^ix, 0) + ^^V(-a:, 0) . (8) 



In a similar manner one may show that 



^{x,3tR/4) 



(1 



^V(x,0) 



(1-z) 



V;(-x,0). (9) 



Thus, a system with an even potential and a purely 
quadratic energy spectrum supports perfect fractional re- 
vivals at quarters of the revival time. Especially interest- 
ing are the fractional revivals at tn/A and Mr/A which 
for an initially localized wave function will consist of two 
perfect, mirrored "cat states". We have failed to find 
perfect fractional revivals at other fractions of the re- 
vival time for this general class of systems (i.e. utilizing 
only parity properties of the eigenfunctions). 

Let us now be more specific, and consider in turn two 
potentials of the type considered above; namely symmet- 
ric cases of the Poschel- Teller (PT) and Rosen-Morse 
(RM) potentials [|l|,0. FT takes the form (with the 
ground state energy set at zero) 



Vsiix) = -A^ + A{A - a)sec^{ax) 



(10) 



defined in the range — 7r/2 < ax < 7r/2, and with energy 
spectrum 



E„^ {A + naf 



A' 



(11) 



In order to have perfect quarterly revivals, it is necessary 
to choose A — Ma, with M an integer. Note that PT 
has an infinite number of bound states, and no scatter- 
ing states. The bound states may be expressed in terms 
of Gegenbauer polynonfials with argument sin(ax). We 
shall restrict our attention to M = 2, in which case the 
energy eigenfunctions are sums of bilinear products of 
elementary trigonometric functions. 



RM takes the form 



Vs2{x) ^ A^ - A{A + a)sech^(aa;) 



(12) 



defined for x on the entire real line, and with energy 
spectrum 



E.„ 



A' 



{A - naf 



(13) 



Again, to ensure perfect quarterly revivals, we choose 
A = Ma, with M a positive integer. RM has only M 
bound states, which may be expressed in terms of the 
Gegenbauer polynomials with argument ismh{ax). 

Although the energy spectra for these potentials are 
not purely quadratic in n, it is a simple matter to redefine 
the quantum number by shifting by Af , in which case the 
perfect quarterly revivals found above take the slightly 
modified form 



1 



1 



^k{x,tR/A) = -(1 -i6i)V'(x,0) + -(1 +z6i)^(-a;,0) 



i^{x,tR/2)^ij{-x,Q) , 
i^{x,itR/A) = \{l + ie)^{x,0) 



(14) 



1 



{i-ie)ilj{-x,Q) 



where 9 = (—1)*^. If one chooses A/ a to be a half- 
integer, the revival at tR/2 is identical to the full revival, 
yet the perfect quarterly revivals are lost. 

Much can be learnt about these systems by preparing 
an initial wave function from a finite number of eigen- 
states, and then studying its evolution using Eq.(g). It 
is common practice [|l| to construct the initial wave func- 
tion using the first N modes, with weights c„ drawn from 
a Gaussian distribution centered at some reasonably en- 
ergetic mode n 



exp 



[n — nf 



2a2 



(15) 



This simulates, for instance, a laser-prepared state in an 
ion trap. We have evolved the wave function in the PT 
and RM potentials using these weights. [We choose equal 
phases for the c„ in order to create a well-localized initial 
wave packet. Choosing random phases does not affect the 
revival structure, but tends to obscure the regularity of 
the quantum carpets.] Aside from confirming the per- 
fect quarterly revivals found above, we have found that 
PT has nearly perfect fractional revival states at ratio- 
nal fractions of the revival time, whereas there is little 
sign of such nearly perfect states for RM (although they 
will appear for much higher energy wave packets ||l[]). In 
Figs.l and 2 we illustrate this by showing the probability 
density p at times t = 0, Ir/^, Ir/A and tfl/3 for PT 
and RM potentials respectively. To test the robustness 
of the quarterly revivals we have also evolved an excited 
wave packet in the RM potential, but away from the ra- 
tionality condition A/a ~ M . We set A/a ~ M + r, 



with r S (0, 1/2). We find imperfect, yet "smooth" frac- 
tional revivals at 1^/2 for general values of r (see Fig. 
3a). However, the fractional revivals at tn/A are more 
sensitive to r, and fade away for r > 0.25 (see Figs. 3b- 
3d). 

The almost perfect fractional revivals for PT may be 
understood intuitively, since for a wave function con- 
structed around a moderately energetic state, the PT 
potential closely resembles the ISW - i.e. the harmonic 
structure of the bottom of the well is barely resolved. 
[One may also make a more quantitative argument by 
changing the basis of Eq. (0) from the PT eigenfunctions 
to the eigenfunctions of the ISW. One finds that the over- 
lap integrals for large n are sharply peaked. Thus the 
dominant part of the wave function may be expanded 
(with the weights c„) in terms of the energy eigenstates 
of the ISW, which, as mentioned before, has perfect frac- 
tional revivals at all rational fractions of tn] Given 
the poor resolution of the structure of the potential, 
one would also expect that the quantum carpets for PT 
closely resemble those found for the ISW. This is indeed 
the case, as shown on the left of Fig. 4, where the charac- 
teristic rays are clearly visible. 



-0.6 -0.4 -0.2 0.2 0.4 0.6 

a) 



■0.6 -0.4 -0.2 0.2 0.4 0.6 

b) ^ 





h . ', . r\ 


- 



-0.2 0.2 0.4 



c) 



■0.6 -0.4 -0.2 0.2 0.4 0.6 

d) '^ 



FIG. 1. Probability density for the PT potential (with pa- 
rameter values a = n, M = 2, N — 30, a — 3.0, and n = 15) 
at times a) f = 0, b) f = ffl,/5, c) t = ifl/4, and d) f = tii/3. 

As an alternative to preparing the initial wave function 
around some energetic state, one can use a weighting 
that favors the lower-lying states (using an exponentially 
decaying distribution for the c„ for example). We have 
studied this case, and indeed the correspondence with the 
ISW disappears, and the evolution of the wave function 
is a slowly modulated classical oscillation (since the well 
has a harmonic minimum) . The almost perfect fractional 
revivals are invisible, and the quantum carpet has no 
structure. The perfect quarterly revivals may still be 
resolved. 
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FIG. 2. Probability density for the RM potential (with pa- 
rameter values a = 1.0, M = N = 20, a = 4.0, and n = 10) 
at times a) t = 0, b) t = tn/b, c) t — ifl/4, and d) f = tn/S. 
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FIG. 3. Probability density for the RM potential (with pa- 
rameter values a = 1.0, M = N = 20, a = 4.0, and n = 10) 
for values of A/ a ^ M + r, with a) r = 0.0 (leftmost), 0.25 
(centre), and 0.5 (rightmost) at time tR/2; b) r = 0.0 at ifl/4; 
c) r = 0.25 at f = tRf-i; and d) r = 0.5 at f = tn/'i. 

As mentioned above, the RM potential shows little sign 
of fractional revivals, apart from the perfect quarterly 
revivals. However, the quantum carpet for this poten- 
tial reveals considerable structure, as shown on the right 
of Fig. 4; the rays of PT are replaced by a complicated 
structure of what appear to be non-linear "world lines" . 
A magnified view of the first perfect quarterly revival is 
shown in Fig. 5. 
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FIG. 4. The quantum carpets (i.e. space-time contour 
plots of the probability density) for PT (left) and RM (right) , 
with parameter values as before, shown for t £ (0,tH/2). The 
darker regions indicate higher probability density. 
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FIG. 5. Detail of the quantum carpet for the RM poten- 
tial around t = tfl/4, showing the emergence and subsequent 
collapse of two perfect cat states. 

We have yet to attain a physical understanding of the 
quantum carpet for RM. Approaches which work well 
for the ISW P] are less useful here due to the compli- 
cated nature of the eigenfunctions. It is interesting to 
note that the structures visible in the ISW/PT quantum 
carpet {i.e. rays) can be interpreted as a superposition 
of coherent wave packets (CWP) which follow classical 
trajectories (with a discrete spectrum of initial velocities 
selected to ensure a perfect revival at t = i/j). Whether 
these CWP can be identified with the coherent states 



p4[ corresponding to the ISW is unclear. The appar- 
ent world-lines of the RM quantum carpet might also be 
interpreted in this way, although the CWP no longer fol- 
low classical trajectories. This is clear from the manner 
in which the world lines proceed through the minimum 
of the potential at i = <fl'/4. 

A computational application of our results is the test- 
ing of numerical algorithms designed to integrate forward 
the time-dependent Schrodinger equation |l5). These al- 
gorithms do not generally rely upon energy eigenfunc- 
tion expansions, and testing them against exact results 
(for general initial conditions, and non-trivial potentials) 
is difficult due to the scarcity of such results in quan- 
tum dynamics. An algorithm which integrates the wave 
function forward in time in the PT or RM potentials, 
and successfully generates {i.e. recovers with good pre- 
cision) perfect quarterly revivals can be trusted in other 
applications. A positive feature of this test is that one 
can implement it for any initial condition (strictly true 
only for PT for which the bound states form a complete 
set). It is interesting to note that algorithms which use 
discrete Fourier modes for free-particle propagation will 
fail to capture perfect revivals, since they are based on 
a Hamiltonian with a discrete lattice Laplacian, thus re- 
placing the pure fc^ spectrum by 2(1 — cosk), although 
they will find increasingly good revival structures as the 
number of Fourier modes is increased. 

In conclusion, we have studied the quantum revivals 
and carpets for systems with a quadratic energy spec- 
trum and an even potential. We have proven that all 
such systems have perfect quarterly revivals, in addition 
to perfect complete revivals. We have studied the time 
evolution of two such systems - the PT and RM poten- 
tials - more closely. From a moderately energetic initial 
distribution of modes, the evolution of the wave func- 
tion in PT (being defined in a finite region of space) has 
many similarities to that of the ISW, with almost per- 
fect fractional revivals at rational fractions of i/?, and a 
quantum carpet with characteristic rays. This similar- 
ity disappears continuously as one decreases the mean 
energy of the wave function, thus allowing better resolu- 
tion of the harmonic minimum of the well. The evolution 
of the wave function in RM shows little fractional revival 
structure, apart from the perfect quarterly revivals. Its 
quantum carpet, although devoid of rays, displays a daz- 
zling pattern, the understanding of which is currently 
being pursued. 

This study has shown that many of the rich dynamical 
properties of the ISW are fairly generic, thus increas- 
ing their experimental relevance. Indeed, the PT and 
RM potentials capture features of real quantum systems, 
which are missing in the ISW. Namely, a spatially varying 
potential energy, with a harmonic minimum; and in the 
case of RM, a finite number of bound states (see also the 
discussion of super-revivals in the finite square well p6|). 
The robustness of the perfect quarterly revivals may well 



be of interest to experimentalists seeking to create per- 
fect cat states from a localized wave packet. Indeed, this 
is the initial entangled state required for quantum com- 
munication |17| (although "entanglement" usually refers 
to two or more degrees of freedom), and which is gener- 
ally created using more complicated laser interferometry. 
(Fabrication of an approximate RM potential may well 
be realisable using semiconductor quantum well technol- 
ogy-) 

Aside from the PT and RM potentials, there are other 
potentials which are isospectral to the ISW (and which 
may be generated using the Darboux transformation 
||13| ). A study of their revival properties may well prove 
worthwhile. As a final remark, it is noteworthy that the 
PT potential (with AI — 2) and the ISW are supersym- 
metric partner potentials |13[| , and therefore share the 
same energy spectrum (bar the lowest state). Whether, 
due to supersymmetry, these systems share other dynam- 
ical equivalents, aside from perfect quarterly revivals, is 
an interesting open question. 
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